N.  M.  Bhatia 

iiid 

Nachbar 

UJ 

ul 


INDENTATION  OF  AN  ELASTIC  MEMBRANE 


BY  A  SPHERICAL  INDENTER 


FEBRUARY 

1964 


PR:rAftEO  FOR  THE  U.S.  ARMY  RESEARCH  OFFICE  (DURHAM) 
UNDER  GRANT  NO.  tJA-AR0|0)3M24  0464 


Department  of  Aeronautics  and  Astronautics 
Stanford  University 
Stanford,  California 


FINITE  INDENTATION  OF  AN  ELASTIC  MEMBRANE 
BY  A  SPHERICAL  INDENTER 

by 

N.  M.  Bhatia 
and 

W.  Nachbar 


SUDAER  No.  182 
February  1964 


Requests  for  additional  copies  by  Agencies  of  the  Department 
of  Defense,  their  contractors,  and  other  Government  agencies 
should  be  directed  to: 

Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22314 

Department  of  Defense  contractors  must  be  established  for  DDC 
services  or  have  their  "need -to -know"  certified  by  the  cognizant 
military  agency  of  their  project  or  contract. 

All  other  persons  and  organizations  should  apply  to  the: 

U.  S.  Department  of  Commerce 
Office  of  Technical  Services 
Washington  25,  D.  C. 


This  work  was  performed  at  Stanford  University 
with  the  sponsorship  of  the 

U.  S.  Army  Research  Office  under  Grant  DA-AR0(D)31-124  G464, 

Project  Number  DA -5 9901007 


ERRATA 


Page 

Line 

Correction 

iv 

4 

For  "Contract"  read  "Contact" 

For  "Indentor"  read  "Indenter" 

1 

22 

For  "is  to  the  membrane  r" 
read  "to  the  membrane  is  r" 

1 

24 

For  "Solution"  read  "A  solution  ..." 

4 

8 

For  "with"  read  "With  ..." 

5 

2 

Should  read:  "equations  of  equilibrium 
[Eqs.  {2.1a,b)  below  with  p  =  O],  ..." 

6 

1 

For  "INDEMTOR"  read  "INDEMTER" 

7 

9 

Should  read:  "The  differential  operator 
transforming  as  [Eqs.  (2.6)  and  (2.7)]" 

7 

10 

Comma  after  equation. 

7 

11 

Should  read:  "Eq.  (2.5b)  can  then  be  written" 

9 

2 

For  "(2.19a,b)"  read  "(2.l4a)" 

Ik 

12 

For  "since"  read  "Since" 

17 

1 

For  "differential"  read  "derivative" 

17 

7 

For  "Eq.  (3.22a)"  read  "(Eq.  3.21a)" 

17 

9 

Should  read  "This,  with  the  use  of  the  left 
Inequality  (3.20),  yields" 

17 

15 

Should  read  "The  total  differential  of  the 
left-hand  side  of  Eq.  (3.12)  yields  an 
expression  for  the" 

18 

14 

For  "C^  ^  Col"  "Co  =  Col" 

21 

8 

For  "yeilds"  read  "yields";  comma  at  end  of  line. 

22 

4 

For  "(3.38a)"  read  "(3.38)" 

22 

17 

For  "P  5  p^"  read  "P  ^  Pl"  . 

Read  this  line  as  "let  Ol  “  ®y  =  1*^^  psi.^ 
and  calculate  =  ;...44,  =  0.0704  lbs,  and 

Fql  =  6.l».  Also  (ciro/®L)®  ' 

=  7-2  X  10'^ 

24 

10 

24 

11 

For  "yeild"  read  "yield"  . 

24 

22 

For  "Limiting"  read  "The  limiting..." 

ABSTRACT 


Rotatlonally  symmetric  stresses  and  deformations  are  considered  for 
a  prestressed  elastic  sheet  of  circular  outer  boundary  loaded  transversely 
by  a  centered  Indenter  with  a  hemispherical  tip.  A  nonlinear  membrane 
solution  is  obtained  for  the  portion  of  the  sheet  that  is  in  frictionless 
contact  with  the  rigid  tip  of  the  indenter.  This  solution  and  the  solu¬ 
tion  that  was  previously  obtained  for  a  prestressed  annular  membrane 
[Ref.  1]  are  used  to  obtain  the  exact  solution  for  stresses  and  de¬ 
flections  in  the  indented  membrane.  Simpler,  limiting  expressions  for 
stresses  and  deflections  are  also  obtained  for  sufficiently  large 
and  sufficiently  small  indenter  loads. 

Comparison  of  computed  results  from  this  analysis  with  the  ex¬ 
perimental  results  of  Jahsman,  Field  and  Holmes  [Ref.  2]  on  stretched 
mylar  membranes  shows  good  agreement  except  in  the  immediate  neighbor¬ 
hood  of  the  indenter.  The  theory  shows  that  yielding  and  plastic 
deformation  of  the  membrane  is  incipient  for  a  relatively  small  value 
of  the  indenter  load,  and,  for  these  experiments,  this  limiting  value 
is  less  than  the  least  value  of  indenter  load  for  which  data  was 
reported . 
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,  combined  loading  parameter  . 
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o  L 


Value  of  F^  for  P  =  P_  . 

Thickness  of  elastic  sheet. 

=  ho  ,  Horizontal  component  of  stress  resultant  at 
ro 

r  =  a . 

Central  indenter  load. 

=  47rEhc  sf  ,  upper  bound  on  P  at  elastic  limit. 
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Indenter . 

Radial  coordinate. 
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Transverse  (vertical)  displacement  component. 

=  -w(o)  ;  Central  deflection. 
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Angle  of  tangent  rotation  • 
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Radial  mid-surface  strain. 
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Poisson's  ratio. 
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Radial  and  circumferential  stress  components. 
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Elastic  proportional  limit  stress. 
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SECTION  1.  INTRODUCTION 


A  theoretical  analysis  has  been  given  by  Nachbar  [Ref.  l]  for  the 
finite,  rotationally  symmetric  deformations  and  stresses  of  a  prestressed, 
annular  elastic  membrane,  or  sheet,  subjected  to  applied  transverse 
loading  only  along  the  inner  edge.  The  sheet  is  initially  flat  and  is 
supported  at  the  outer  edge  where  stretching  is  also  applied.  An  ex¬ 
plicit  solution  was  obtained  for  the  case  of  transverse  loading  introduced 
through  a  rigid  plug  or  disk  which  is  attached  centrally  to  the  annular 
membrane.  Experimental  results  for  a  similar  problem  have  also  been 
given  by  Jahsman,  Field  and  Holmes  [Ref.  2].  However,  the  experiments 
used  a  complete  membrane  ana  a  rigid  Indenter  with  a  hemispherical  tip 
in  order  to  apply  transverse  load.  Thus,  stresses  and  deformations  of 
the  membrane  in  the  immediate  neighborhood  of  the  indenter  tip  are  to 
be  expected  to  be  different  from  those  predicted  by  the  plug  analysis. 

In  the  present  paper  is  considered  the  problem  of  finite,  rotationally 

symmetric  deformations  and  stresses  in  a  prestressed  circular  flat  sheet, 

of  outer  radius  a  ,  due  to  transverse  loading  at  the  center  by  a  rigid 

indenter  with  a  hemispherical  tip  of  radius  c  .  This  problem  will 

henceforth  be  referred  to  as  the  Indenter  problem.  A  small  strain, 

elastic  solution  for  the  Indenter  problem  is  obtained.  Figure  1  shows 

the  indenter  problem  geometry  and  the  nomenclature.  The  radial  distance 

is  to  the  membrane  r,  and  r  =  b  denotes  the  point  of  tangency  of  the 

sheet  with  the  indenter:  b  is  a  function  both  of  load  P  and  of  H  . 

o 

Solution  for  the  portion  of  the  sheet  (O  5  r  S  b)  in  frictionless  con¬ 
tact  with  the  indenter  is  obtained  in  Section  2.  This  portion  will  be 
called  the  constrained  region.  This  solution  and  the  results  for  the 
free  region  (annular  portion  of  the  membrane  which  is  not  in  contact 
with  the  Indenter)  are  used  in  Section  3  to  obtain  the  complete  solution 
to  the  Indenter  problem. 

The  results  for  the  free  annular  region  [Ref.  1]  which  are  essential 
for  the  present  paper  are  summarized  below.  These  are  valid  for  small 
strain  and  moderate  rotation  restrictions.  The  inner  radius  is  b  and 
the  outer  radius  is  a.  Displacements  normal  and  parallel  to  the  initial 
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RIGID  INDENTER  WITH 
HEMISPHERICAL  TIP 


I  FREE  REGION  (bsrSa) 

II  CONSTRAINED  REGION  (OSrSb) 


DEFORMED  EQUILIBRIUM  CONFIGURATION 

Fig.  1.  Deformation  of  the  Elastic  Sheet  —  Deformed  Equilibrium 
Configuration. 


t. 


DEFORMATION  OF  A  MEMBRANE  ELEMENT  IN  THE 

CONSTRAINED  REGION 


Fig.  2.  Deformation  of  the  Elastic  Sheet  —  For  the  Constrained 
Region. 


3 


plane  are  w  and  u  ,  respectively.  The  angle  of  rotation  of  the 

tangent  to  the  midsurface  from  the  initial  plane  is  P  .  The  sheet  has 

uniform  thickness  h  and  elastic  coefficients  E,V  .  Nondimensional 

stress  variables  s  and  s.  are  defined  as 

r  0 


and  the  independent  variable  is  chosen  as 


with  use  of  the  nondimensional  shape  parameter  €  and  the  load 
parameter  p  , 


^  b  ^  P 

^  ”  a  ’  ^  "  27rEhb 


(1.1c) 


a  nondimensional  stress  function  F  is  defined  as 

F  =  (Uep)^/^  ctn  p  s  (4€p)^/^  y^/^  (l.ld) 

Stresses  and  displacements  are  given  in  terms  of  F(y)  by  Eqs.  (3.7)» 
(3.8),  (3.9),  (3.11)  and  (3.13)  of  Ref.  1  as  follows: 


s 

r 


w 

a 


P  - 


(4€p) 


1/3 


dy 


(1.2) 

(1.3) 

(1.1) 

(1.5) 


-  1  - 


(1.6) 


Equations  (1.3)  and  (1.4)  are  seen  to  satisfy  identically  the  differential 
equations  of  equilibrium  (Eqs.  2.1a,b  below  with  p  =  O),  and  the  equa¬ 
tion  of  compatibility  [Eq.  (2.5b)]  becomes,  with  the  use  of  Eq.  (1.6), 
a  differential  equation  for  F(y): 

^  +  2F"^  =  0  (1.7a) 

dy^ 


The  integral  of  this  corresponding  to  s-  >  0  has  the  general  form 


F(l) 

F(y) 


(1/2) 


dv  , 


(1.7b) 


where  C  is  a  real-valued  constant. 

0 

The  boundary  condition  at  y  =  1  is  Ehs  =  ho  =  H  .  Thus  H 

r  ro  o  o 

determines  the  prestress  .  The  parameter  F^  expresses  the 

combined  loading  condition,  through  use  of  Eq.  (1.4),  as 

2/3 

2 

The  boundary  condition  at  y  -  €  is  expressed  as 

F(€^)  =  F^  (1.8b) 

where  F^  is  to  be  determined  from  conditions  of  continuity  with  the 
constrained  region.  When  €,F  ,F  are  given  so  that  F  >  F  >  0  ,  C 
is  uniquely  determined  from  Eq.  (l.7b)  expressed  at  y  =  as  follows: 


(1.9) 
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SECTION  2.  SHEET  IN  FRICTIONLESS  CONTACT  WITH  INDENTOR 


Finite  deformation,  small  strain  equations  for  the  constrained 
region  (0  S  r  £  b)  of  the  sheet  are  [Ref.  3]  those  of  equilibrium, 


d 

(rs  cos  S)  -  8^  +  ^  sin  3  =  0 
dr  '  r  0  Eh 


(2.1a) 


^  (rs  sin  3)  -  cos  3=0 
dr  '  r  ^  Eh 

elasticity  [0  <  v  £  (l/2)]  , 


=  s  -  Vs 
r  r  o 


(2.1b} 

(2.2a) 

(2.2b) 


and  strain-displacement, 

(l  +  €  )  sin  P  =  ^ 
r  dr 


(l+€  )cos3=l+f^ 

'  r  dr 


(..^3a) 

(2.3b) 

(2.3c) 


When  p  is  eliminated  between  Eqs.  (2.1a,b),  there  obtains 


d?  -  ^9  ® 


(2.M 


The  compatibility  relation  obtained  by  combining  Eqs.  (2.3b,c), 


—  (rCg)  +  1  -  (1  +  e^)  cos  3  =  0  , 


(2.5a) 


can  be  written  in  a  more  convenient  form,  with  the  use  of  Eqs.  (2.2) 
and  (2.U),  as 


r 


dr 


+  s  )  +  (1  +  s  ^  s  )(l  -  cos  3)  =  0 
r  or 


(2.5b) 
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Since  the  membrane  is  assumed  to  deform  smoothly  onto  the  indenter,  the 
geometric  relation  for  the  constrained  region 


sin  0  =  =  -  (1  +  €  )  (2.6) 

c  c  '  6 

is  apparent  from  Fig.  2. 

It  is  next  shown  that  f  0  |  must  necessarily  be  of  the  order  of  the 
square  root  of  the  maximum  elastic  strain  or  smaller.  From  Eq.  (2.3b) 
and  the  derivative  of  both  sides  of  Eq.  (2.6), 


d0 

dr 


-1 

c 


(2.7) 


transforming  as  Eqs.  (2.6)  and  (2.7)*  The  differential  operator 


r  iQ  =  I 

dr  c  d(r/c )  d0  \  ^  ^0 


sin  0 


ill . 

dg 


Equation  (2.5b)  can  be  written 

(l  +  sin  0 

Because  of  the  assumption  of  small  strains,  js^|  «  1  and  |Sq|  «  1  , 
the  square  bracketed  terms  in  the  above  equation  can  be  approximated  by 
unity,  and  the  equation  becomes 


d0 


+  tan  I 


0 


(2.8) 


Let  s(0)  =  s^  +  Sq  ;  then  the  integral  of  Eq.  (2.8)  can  be  written  as 


s(0) 


s(0)  +  2  in  cos 


0 

2 


(2.9) 


or  also  as 

cos  I  =  exp  i  [s(0)  -  c(0)]} 


(2.10) 


Equation  (2.10)  shows  that  s(0)  ^  s(0)  for  all  0  >  0  .  Furthermore, 
since  |s|  is  bounded  by  some  linear  elastic  limit  value  denoted  by 

-  7  - 


(2s  ),  where  s  «  1  ,  then  Eq.  (2.10)  implies 
L  L 

[1  -  cos  (3/2)]  <  1  -  exp(-s^)  i  Sj^ 

2  ***  o 

so  that  3  S  os  .  Therefore,  if  s  is  sufficiently  small,  it  is 

L  L 

necessary  to  assume  as  consistent  with  the  small  strain  approximation 
that  the  rotation  is  finite  but  "moderate",  viz. 


2 

3  «  1 

Hence  Eqs.  (2.6)  and  (2  9)  are  approximated  as 


(2.11) 


sin  3=3=— 
c 


(2.12a) 


s(3)  =  s(0)  -  ^  3^  (2.12b) 

If  Eq.  (2.4)  is  written  in  the  form 
ds 

r  +  2s  =  s  -  S-(l  -  cos  3) 
dr  r  0 

then  with  the  use  of  Eqs.  (2.12a,b)  and  the  small  strain  assumption, 
this  becomes 
ds 


P  ^  .  23_.  .  s{0)  -  1  £>2 


(2.13) 


The  solution  is 


.  4 . 4^  .  i,  e2 


r  p2  ■  2  -  15 


(2.lUa) 


and  from  Eq.  (2.12b)  , 

^2  siq)  5  _2 

^0  =  -  ^  ^  2  -  ^ 


(2.l4b) 


With  use  of  these  relations  in  Eqs.  (2.2b)  and  (2.3c),  u  is  determined 
to  be 

C, 


^  =  -(1  4  v)  ^  +  I  (1  -  v)3  s(0)  -  ^  (3  -  v)3^ 


(2.15) 
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The  contact  pressure  p  between  the  membrane  and  the  indenter  is 
determined  by  using  Eqs.  (2.1b),  (2.12a)  and  (2.19a,b): 


Eh 

P  =  - 


s(o)  -  i 


(2.16) 


An  expression  can  be  given  for  s(0)  in  terms  of  >  where 

=  (b/c)  ,  by  using  Eq.  (2.lUa): 


6(0)  =  2  6^{Pj^)  .  2  ^  1  p2 


{2-i7) 


Equations  (2. l4 )  to  (2.1?)  are  valid  for  ^  0  ,  but  if  the  membrane 

is  to  be  elastic  for  all  3^0,  then  =  0  .  For  the  case 

=  0  ,  Eqs.  (2.IU)  to  (2.16)  become  with  use  of  Eq.  (2.17): 


(2. 18a) 


(2.18b) 


ul£l  ^ 


3 


Eh 


(1  -  V)  el  .  pS 


p{e)  -  f 


1  .2  1  „2 


'•'  ^  ^  ^  -  IT  ^ 


(2.18c) 

(2.l8d) 


The  displacement  w(0)  is  found  by  integration  of  Eq.  (2.3a): 


w(3)  =  w(3^)  -  I  (3^  -  3^) 


(2.l8e) 


A  necessary  condition  to  be  satisfied  for  contact  between  indenter 
and  membrane  in  0  5  3  S  ,  is  that  p(3)  ^  0.  This  coiidition  is 
seen  to  be  fulfilled  if 


SgOj,)  >  0 


(2.19) 


This  latter  condition  is  necessary,  however,  for  the  stability  of  the 
free  membrane  against  wrinkling  [Ref.  1,  Eq.  (3.l6)]  and  must  be  imposed 
for  a  rotationally  symmetric  solution  to  the  indenter  problem.  It  will 
be  shown  in  the  following  section  (see  Eq.  (3.39)]  that  (2.I9)  is 
satisfied . 


_  o  _ 


The  equilibrium  of  the  indenter  at  the  point  of  tangency,  0 
requires 


P  =  27rEh  bp,  S  (P,  ) 
b  r  b 


which  is  written  as 


P 

27ifEha 


(2.20) 


iO 


SECTION  3.  SOLUTION  OF  THE  INDENTER  PROBLEM 


Solutions  from  Section  1  for  the  free  region  and  from  Section  2 

for  the  constrained  region  are  now  used.  Solutions  for  these  two  regions 

are  distinguished  by  the  use  of  y  as  independent  variable  in  the  free 

region  and  of  S  as  independent  variable  in  the  constrained  region.  At 

2 

the  point  of  tangency  r  =  b  (y  =  €  for  the  free  region  and  S  =  3^ 
for  the  constrained  region),  conditions  of  continuity  of  displacements 
and  normal  stress  are  imposed.  From  Eq.  (l.^), 


where,  from  Eq.  (l.lc). 


P 

"  27rEha 


(3.2) 


Continuity  of  s^  requires 


s 

r 


)  =  s^(P,) 


(3.3) 


This  equation,  after  substitution  from  Eqs.  (1.8b),  (2.20),  (3.1)  and 
(3.2),  can  be  written  as 


F  = 

€ 


TrEha  / 


c 

a 


(3.4) 


Note  that  F^  does  not  depend  upon  e  . 

Continuity  of  the  displacement  u  is  invoked  to  determine  b,  or 
what  is  equivalent,  the  nondimens ional  parameter  e: 


u(€^)  =  u(3^) 


Equations  (1.2),  (l.7a),  (l.8b)  and  (3.2)  give 


u(c^) 

b 


(3.5) 


(3.6) 


Equation  (2.l8c)  gives 


11  - 


(3.7a) 


=  (1  -  V)  S^(p^) 


2 

b 


which  can  be  written  in  the  convenient  form,  using  Eqs.  (3.1)  to  (3.4), 


u(B^) 


P 

47rEha  / 


(3.7b) 


Substitution  from  Eqs.  (3.7b)  and  (3.6)  into  Eq.  (3.5)  leads  to  the 

2 

following  equation  quadratic  in  e  , 


4 

€ 


2 

€ 


0 


(3.8) 


The  only  positive  root  is 


2 

c 


r  f 


4?" 


1/2 


which  can  be  written  simply  as 


€ 


(3.9) 


Equations  (1.9)  and  (3.9)  are  used  to  determine  €  and  C  with 

o 

and  ,  or  equivalently,  with  P  and  ,  through  the  use  of 

Eqs.  (l.8a)  and  (3.4),  considered  as  independent  variables. 

Certain  characteristics  of  the  solutions  e  and  to  Eqs.  (1.9) 

and  (3.9)  for  fixed  H  >  0  are  now  studied.  Restrictions  of  these 

o 

developments  to  small  indenter  radii  will  be  made  subsequently  as 
sufficient  for  applied  purposes  and  in  the  interest  of  brevity. 


F 

o 


First,  it  is  observed  that  F^  is  expressible  as  a  function  of 

and  H  by  means  of  Eqs.  (l.8a)  and  (3.9) 
o 


Therefore,  for  fixed  H  ,  F  may  be  regarded  as  a  function  of  F 

o  €  o 

only,  viz.  •  Substitution  of  into  Eq.  (3.9)  yields 
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the  function  €  =  .  Then  the  right-hand  side  of  Eq.  (1.9)  can 

be  defined  by  the  function  l(F  ,C  )  as 

'  o  o 


o  0 


(3.11) 


and  then  Eq.  (1.9)  can  be  written  as 


€^(F  ,C  )  +  I(F  ,C  )  -  1  =  0  (3.12) 

''  0  o  o  o 


Equation  (3.12)  will  be  shown  to  determine  implicitly  a  function 

C^(F^^)  I  which  can  be  regarded  as  the  solution  to  the  problem.  Since 

0  ^  €  <  1  is  a  geometrical  constraint,  then  by  Eq.  (3.12) 

0  <  I[F  ,C  (F  )]  i  1.  It  is  necessary  to  consider  the  left-hand  side 
0  0  0 

of  Eq.  (3.12)  in  the  plane  within  the  allowable  domain  which 

is  determined  by  the  two  inequalities 


+  C  >  0 
0 


0  <  F  < 
€ 


1/3 


2/3 

<  F 

o 


(3.13a) 


(3.13b) 


Inequality  (3.13b)  follows  from  Eq.  (3.10)  and  the  inequality  I  >  0, 
which  implies  ^  ^  ‘ 

It  is  evident  from  Eqs.  (3.9),  (3.10)  and  (3.1l)  that 

and  I(F  ,C  )  are  continuous  functions  of  F  and  C  in  the  allowable 
0  0  0  0 

domain.  The  partial  derivative  with  respect  to  of  the  left-hand 

side  of  Eq.  (3.12)  is 


*) 

It  is  shown  in  Ref.  1  that  if  wrinkling  is  to  be  avoided, 
Co(f’o)  roust  be  confined  to  a  smaller  domain  which  lies  wholly  within 
the  allowable  domain. 
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(3.lUa) 


-(1/2),^  x-(l/2) 


-3/2 

-/ 


Therefore,  in  the  allowable  domain,  the  right-hand  side  above  is 
negative,  and  so 


3r  *  I  - 1)  <0 


(3.l4b) 


This  partial  derivative  is  a  continuous  function  of  F  and  C  . 

o  o 

At  this  point,  the  further  restriction  to  small  indenter  radii  is 
made,  viz. 


-  <  <  1 
a 


(3.15a) 


From  this  it  necessarily  follows  that 
€  <  <  1 

since  a  ^  o.  and 
ro  L 


^  ^  <<  I 

then  E  ^  o  <  <  1  ,  and  therefore 
ro 

(H  /Eh)  <  <  1 
o 


(3.15b) 


(3.15c) 


(3.15d) 


Consequently,  it  can  easily  be  shown  from  Eq.  (3.10)  and  inequalities 


(3.13b),  (3.15a,d)  that 


(3.15e) 


Ik  - 


When  C  =  0  ,  let  F  =  F*  be  the  root  of  Eq,  (3.12).  Tnen  with 
o  o  o 

the  use  of  Eqs.  (3.9)  and  (3.1l),  Eq.  (3.12)  yields 


„«(3/2)  _  3 
o  2 


(3.l6a) 


In  view  of  (3. 15e) 


o 


(3.16b) 


It  follows  from  inequalities  (3.13),  (3.15a,d)  that  the  point  F^  =  F* 

C  =  0  ,  lies  within  the  allowable  domain, 
o 

With  these  restrictions,  we  will  now  prove  the  existence  of  C 

o 

as  a  continuous  function  of  F^  in  the  following  subdomain  of  the 
allowable  domain: 


C  2  0  and  F  §  F*  (3.1?) 

2 

Inequality  (3.l4b)  Implies  that  (e  +  I  -  l)  is  a  monotone,  strictly 

decreasing  function  of  C  for  a  fixed  F  .  Also  for  C  =0  and 
°  o  o  o 

F  >  F"*  ,  it  can  readily  be  shown  that 
o  o 

+  I  -  1  >  0  (3.18a) 


It  can  also  be  shown,  with  ■'he  use  of  Eqs.  (3.9)  and  (3.11),  that 


lim 

Co'-® 


F  fixed 
o 


2 

€ 


0 


and 


lim  1=0  , 

Co-® 

F  fixed 


from  which  it  follows  that 


lim  (e^  +  I  .  i)  =  .1 

Co-  ® 

F  fixed 
o 


(3.18b) 


Therefore,  for  each  ^o  ^  ^o  '  follows  from  (3.l4b),  (3.l8a.b)  that 
there  exists  a  unique  positive  value  of  C^  ,  such  that  Eq.  (3.12)  is 
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satisfied.  The  continuity  of  C^(F^)  follows  from  a  theorem  on  implicit 
functions  (Ref.  4). 


In  the  subdomain  (3.17)>  upper  and  lower  bounds  for 


obtained  when  Eq.  (3.1l)  is  rewritten  as, 

r 


I(F  ,C  )  = 
o  0 


r 


dv  - 

o  o 


F  (F  ) 

e  o 


0 

r 

J 

1  - 

i  -1 

dv  > 


V. 


J 


This  becomes 

)  = 


P 

O'  O'  ^1/2  o 


(3.19a) 


Where 


,-(1/2) 


1  -  1  t 


vC. 


dv 


(2.19b) 


It  is  easily  shown  that 

-1 


0  <  R  <  C 


,£n 


1  +  C  F 
_ o  o 

1  +  C  F 

O  € 


(3.19c) 


Then  Eqs.  (3.12)  and  (3.l6a)  are  used  to  obtain 


( 1  -  =  F  -  F  -  R  <  F  , 

O  O  €  O 


v/hich  yields 


'  -  - 

O 


€  1 


F  C 

O  O 


lr\ 


1  +  F  C 
o  o 

1  +  F  C 

€  O 


<  C  <  F^[l  +  O(e^)  ]  (3.20a) 


o  o 


With  the  right  Inequality  and  Eq.  (3.10),  an  upper  bound  for  is 

obtained 

1/2 


F  C  <  2(  -  F^/^  [1  +  0(€^)  ] 

€0  \  ^  /  a  o 


-  l6  - 


(3.20b) 


The  partial  differential  of  the  expression  in  Eq.  (3.12)  with 


respect  to  is 


=  (f 

o  \  o 


,-(l/2) 


*  3r 

o 


F  ^  F  ^o 
€  \  € 


L 


(3.21a) 


where  with  the  use  of  Eq.  (3.10), 


hF  -  . 

r— ^  =  -  -  F  F"^ 

5f^  2  €  o 

o 


(3.21b) 


Then  Eq.  (3.22a)  becomes 


^  (€^  +  I  -  1)  =  /i-  +  C 


sr 

o 


F  o 
o 


•(1/2)  £  \  (2  ^ 
■  Fo  "  ^o  (^€ 


(1/^) 


(3.21c) 


This  with  the  use  of  the  left  inequality  (3.20)  yields 


r 


(^2  .  I  -  1)  >  f  ^ 

O  ''  ' 


l+0(€^)+i- 

F"^ 

O 


-d/fe)  ,  v-{l/2) 


V 


(3.22a) 


From  this  Inequality  it  is  evident  that 


^  (e^  +  I  -  1)  >  0 


W 

o 


(3.22b) 


in  the  subdomain  (3.17). 

The  differential  of  Eq.  (3.]2)  yields  an  expression  for  the 

derivative  of  C  (f  ) : 

o  o 


dC 

c 

dF 


o 


^  (c^  +  I  -  1) 


o 


3r  ^ 


(3.23) 
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In  view  of  inequalities  (3.l4b)  and  (3.22b),  this  derivative  is  positive 

for  all  F  >  F*  .  This  implies  that  the  solution  to  Eq.  (3.12)  ex- 
0  0 

pressed  by  C^(F^)  Is  a  continuous,  monotone  strictly  Increasing  function. 

Approximate  solutions  for  the  Indenter  problem  will  next  be  obtained 
for  certain  limiting  cases.  Consider  first  the  case  when 


P 


^TrEhc  s 


2 

L 


and 


Then,  with  the  use  of  Eq.  (l.8a) 


F 

o 


P=P 


L 


(3.25a) 


(3.25b) 


(3.25c) 


It  can  be  seen  with  the  use  of  inequality  (3.25b),  that  F  >  F*  ,  and 

O  Jj  o 

therefore  the  solution  for  P  =  P^  is  in  the  subdomain  (3.17).  Upon 

L 

substitution  of  F^  =  F^^^  into  the  right-hand  inequality  (3.20b),  there 

is  obtained  for  F  =  F  ^  and  C  *=  , 

€  €L  O  OL 


F  C  -  F  ,  <  2 
e  o  €L  oL 


ro 


[1  -  0(€‘^)] 


(3.25d) 


Therefore,  if  the  ratio  (o  /o^)  obeys 

ro  L 


2 


ro 


«  1 


o 


(3.26a) 


where  the  left  Inequality  is  obtained  from  Eq.  (3.25b),  then  it  follows 

that  for  F  =  F  . 

o  oL 


F  C  =  F  <<  1 

€0  €L  oL 


(3.26b) 


Hence  the  following  order  relation  will  hold,  by  continuity,  for  F^ 

in  some  neighborhood  of  F  and  for  P  in  some  neighborhood  of  P  : 

OL«  ]_< 


-  18  - 


0  <  F  C  «  1 
e  o 


(3.26c) 


Let  Eq.  (3.9)  be  approximated  for  amall  as 


^  2i  -  1)^^^ 


[1  -  O(F^C^)] 

£  I  P 

a  \  2T7Eha  / 


(3.27a) 

(3.27b) 


Hence 

b  a  2(V^ 


(3.27c) 


Expressions  for  3^  ,  s^(3^)  and  s^(0)  for  P  near  are  obtained 

from  Eqs.  (2.l8a)  and  (2.19)  with  the  use  of  the  approximation  (3.27c): 


1/4 


%(V2  .  1)  V2'®'C 


) 


1/2 


l6{^/2  -  1) 


3 


(3.28a) 

(3.28b) 


s  (0) 
r 


(3.28c) 


It  is  noted  that  these  relations  are  independent  of  the  prestress  and 

membrane  radius  a.  Also  the  maximum  stress  in  the  membrane,  which 

occurs  at  P  -  0(i.e.  r  =  O),  increases  with  P.  Then  for  P  =  P.  ,  the 

Li 

maximum  stress  in  the  membrane,  evaluated  with  the  use  of  Eqs.  (3.28a,c) 
and  (3.25a),  is 


s  (0) 
r 


(3.29) 


which  is  the  clastic  limit  or  yield  stress  0^  .  Therefore,  since  the 
stresses  are  limited  by  a  ,  it  is  concluded  that  if  the  prestress  is 

Ij 

within  the  limit  given  by  (3.26a),  then  P  represents  an  upper  bound 

L 

on  P  ,  viz. 

P  5  P^  ,  (3.30) 
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If  P  is  sufficiently  small,  so  that 


F  »  1  , 
o 


(3.31a) 


then  it  is  easily  shown  with  the  use  of  inequality  (3.20a)  that 


C  « 

o  o 


(3.31b) 


Sufficiently  large  values  of  F^  are  now  considered,  such  that 


F  C  »  1  and  F  F^  »  1  . 

€  O  €  O 


(3.31c) 


This  inequality,  with  the  use  of  Eqs.  (l.8a)  and  (3.4),  implies  that 


E _  P 

(H  /hf 
o 


«  1 


(3.31d) 


With  the  use  of  iiiequality  (3.31c),  Eq.  (3.9)  can  be  approximated  as 


^  ‘""c  r  3 


(3.32a) 


This  expression  with  the  use  of  Eqs.  (3.31b),  (l.8a)  and  (3.10)  yields 


c  /  P 


o  2  I  H 
27ra  \  o, 


(3.32b) 


An  approximation  for  the  integral  I(F^,C^)  in  Eq.  (3.1l)  is  also 
obtained  with  the  use  of  Inequalities  (3.31a,c)  as 


(^  ■  r  •  5fV  r 

C  '  \  o  O  O  € 

o  ' 


(3.33) 


An  expression  for  deflection  w(€  )  is  obtained  from  Eqs.  (1.5)  and 
(1.7a); 


,  2v  /H  /Eh 

w.(5.  .)  _  .2  _2l_ 
r  If 

\  o 


F  F-  +  C 
o^Fo  o 


F  ^  +  C 

€\F^  O 


Cod  -  c  ) 


(3.34a) 
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with  the  use  of  Eqs.  (3.32p'  and  (3.33)  and  inequalities  (3.31a,c), 
this  equation  is  approximated  as 


/h  \ 

W(€  )  ^  O 

a  “  "\Eh / 


H 


(PC  )'^ 


(3.34b) 


o  o 


With  the  use  of  Eqs.  (3.31b),  (l.8a)  and  (3.10),  Eq.  (3.34b)  yields 

=  -  ilir  [¥^  (r)] 

The  central  deflection  w(3  =  O)  is  obtained  from  Eq.  (2.l8e)  as 


w(0)  =  .  I  ^  w(€^) 


(3.35a) 


This  yellds,  with  the  use  of  Eq.  (3.34c)  and  (3.32b) 


w(0)  =  - 


1  +  in 


2^2  "o' 


(3.35b) 


Expressions  for  obtained  with  the  use  of 

approximation  (3.32b)  and  Eqs.  (2.l8a)  and  (2.20)  : 

s_.(ej,)  -  ^  (3.J6a) 

H 


®r^°^  Eh  327rc  H 

o 


(3.36b) 


The  final  topic  considered  in  this  section  is  the  stability  of  the 
free  membrane  against  wrinkling.  It  has  been  shown  in  Ref.  1  (see 
Appendix;  II,  p.  43),  that  the  wrinkling  stability  condition  is  satisfied  if 

60(1)  >0  (3.37a) 

and 

S0(€^)  >  0  (3.37b) 

It  can  easily  be  shown  that  (3.37a)  is  satisfied  if 
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(3.38) 


C 

o 


1 


-r  T 

o 


This  inequality  is  satisfied  both  for  =  F*  and  for  F^  »  1  . 

Furthermore  the  numerical  solutions  for  C  (F  )  show  that  inequality 

o  o 

(3.38a)  is  satisfied  for  all  F  >  F*  . 

'  o  o 

2 

An  expression  for  s-(€  )  is  obtained  with  the  use  of  Eqs.  (2.l8b) 
and  (2.20)  as 


kP  I 

TTEhc  \  a  j  ^ 


(3.39) 


With  the  use  of  Eq. 

and  2  <  F  C  <  “  , 
e  o 


(3.0),  it  can  be  shown  that  for  both 


2 

€ 


< 


2(^2'.  1)F^''^ 


0  <  F  C  <  1 

€  O 


(3.i;0a) 


Furthermore  from  numerical  solutions  this  inequality  has  been  observed 
to  be  satisfied  for  1  <  <  2  .  Inequality  (3.40a),  with  the  use  of 

Eq.  (3.4),  implies  that 


< 


(3.40b) 


Therefore,  with 
that  inequality 
particular,  the 


the  use  of  Eqs.  (3.39)  and  (3.40b),  it  is  concluded 

(3.37b)  is  satisfied  for  ail  F  >  F*  •  Hence,  in 

o  o 

membrane  is  stable  against  wrinkles  for  P  5  . 


SECTION  h. 


NUMERICAL  RESULTS  AND  COMPARISON  WITH  EXPERIMENTS 


Experimental  results  have  been  reported  (Ref.  2)  for  the  central 

load-transverse  deflection  characteristics  and  the  deflection  profile  of 

*) 

a  mylar  sheet  stretched  in  its  plane  by  dead  weight  loading.  This 

load  is  referred  to  as  the  platen  load.  The  sheet  was  supported  by  a 
ring  of  5-ln  inner  radius.  The  transverse  load  P  was  applied  at  the 
center  of  the  sheet  by  a  load  probe,  or  indenter,  having  a  hemispherical 
tip  of  l/l6-in  radius.  The  membrane  deflections  were  measured  both  by 
a  dial  gage  and,  in  the  immediate  vicinity  of  the  Indenter,  by  measure¬ 
ments  from  photographic  enlargements.  Stresses  were  not  measured  directly, 
but  o^  was  calculated  from  an  equilibrium  equation  [the  present  Eq. 

(2.1b)  with  p  =  0],  and  values  of  S  were  obtained  by  numerical  differ¬ 
entiation  of  the  deflections. 

The  procedure  used  to  compute  numerical  results  from  the  present 
theory  is  first  to  obtain,  for  given  values  of  P  and  (l.e.,  for 

a  given  Fq),  the  root  C  of  Eq.  (3.12).  Since  for  these  experiments, 
(c/a)  =  0.0125,  inequalities  (3.15a,b)  follow.  Tlien,  for  each  >  F* 

[Eq.  (3.16)],  the  root  ^  ^  obtained  using  numerical  integration 

to  evaluate  l(F^,C^),  Eq.  (3.il).  Then  c  -  6(F^,C^),  follows  immediately. 

With  the  use  of  these  values  of  e  and  C  ,  and  with  application  of 

Eqs.  (1.7b),  (l.i+),  (1.3),  (1.2)  and  (i.5),  stresses 

flections  w,u  in  the  free  region  of  the  membrane  are  calculated. 

Stresses  and  deflections  in  the  constrained  region  of  the  membrane  are 
then  determined  with  the  use  of  Eqs.  (3.3),  (2.18)  and  (2.20);  note  that 

p 

=  (a/c)€  ,  and  that  w(€  )  from  Eq.  (1.5)  is  equal  to 
Eq.  (2.l8e). 

The  values  of  0^^  could  not  be  obtained  experimentally  to  within 
less  than  15^  error  (see  Ref.  2,  Table  Al),  probably  because  of  friction 
between  the  sheet  and  the  outer  supports.  Even  for  constant  platen 


*) Properties  of  mylar  sheet;  Thickness  h  -  6.0  x  10^  psi 
Poisson's  ratio  v  =  0.3  ,  Yield  stress  in  uniaxial  tension 
0  =  10^  psi. 

*7 
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is  a  weak  function  of  P.  To  compare 


loads  in  the  experiment, 

predictions  of  the  present  theory  with  the  experiments,  the  numerical 

value  of  o  for  each  value  of  P  is  determined  by  adjusting  o 

ro  ro 

in  the  theory  so  that  the  predicted  deflections  near  the  outer  edge 
r  =  a  are  in  agreement  with  the  experimental  deflections.  The  values  of 
thus  found  agree  roughly  with  the  experimental  estimates. 

In  Figs.  3  and  k  are  shown  curves  for  the  principal  stresses  and 

the  transverse  deflection  w  as  predicted  by  theory  for  P  =  0.66  lbs 

(300  grams)  and  a  calculated  value  of  =  fi50  psi;  for  these  values, 

F  =  1.4U,  P-  =  0.0704  lbs,  and  F  .  =6.4.  Also  (0  /o_ =  7.2  x  lO” 

o  L  oL  ro  L 

hence  P  =  P.  at  the  yeild  limit,  in  view  of  inequality  (3.25c)  and 

Li 

Eq.  (3.29).  Therefore,  since  P  >  P.  and  F  <  F  .  ,  the  stresses  in 
the  neighborhood  of  the  indenter  as  predicted  by  the  theory,  exceed  the 
yield  limit  of  the  material  (See  Figs.  3  and  4).  This  indicates  that, 
plastic  deformation  of  the  membrane  in  the  neighborhood  of  the  indenter 
must  have  occurred  in  the  experiment  causing  larger  transverse  displace¬ 
ments  w  than  those  predicted  by  the  theory. 

In  Fig.  3>  central  load-def lectior  characteristics  are  shown  in 
curves  1,  2,  3  —  each  for  a  fixed  prestress.  The  values  of  0^^ 

used  are  those  determined  as  mentioned  above  for  platen  loads  of  10,20 
and  40  lbs.  with  the  respective  Indenter  loads  of  300,  350  and  400  gms. 
Limiting  value  of  load  P^  =  0.0704  lb.  is  also  indicated.  Since  ex- 

ij 

perimental  data  was  not  available  for  loads  below  p  ,  no  direct 

L 

comparison  between  theory  and  experiments  was  possible.  However,  dis¬ 
crepancy  between  theory  and  experiment  decreases  for  smaller  loads. 

In  the  neighborhood  of  limit  load  P  ,  theory  and  experiments  show 

L 

good  agreement. 


LATERAL  DEFLECTION  -w (inches)  STRESSES  o*x|0  (psi) 


40 


PROBE  LOAD  P«0.66  lbs 
0-0  EXPERIMENTAL  DATA 


Fig.  3.  Deflection  w  and  Stresses  a  and  vs.  (r/a)  for  the 

Membrane . 


(psi) 


r  (inches) 


Fig.  4.  Deflection  w  and  Stresses  and  vs.  (r/a)  for  the 

Membrane  in  the  Neighborhood  of  the  Indenter. 
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o  o  o  EXPERIMENTAL  DATA 
-  COMPUTATION  FROM  THEORY 


O 


w 
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INDENIER  LOAD  P(lbs) 

Fig.  5-  Central  Deflection  vs.  Indenter  Load  P  for  Fixed 


REFERENCES 


1.  W.  Nachbar,  "Finite  deformation  of  a  prestressed  elastic  membrane", 
SUDAER  No.  lUl,  November  I962. 

2.  W.  E.  Jahsman,  F.  A.  Field  and  A.  M.  C.  Holmes,  "Finite  deformations 
in  a  prestressed,  centrally  loaded,  circular  elastic  membrane". 
Proceedings  of  the  fourth  U.S.  National  Congress  of  Applied 
Mechanics,  ASME,  I962. 

3.  E.  Relssner,  "Rotatlonally  symmetric  problems  of  thin  shells  of 
revolution".  Proceedings  of  the  third  U.S.  National  Congress  of 
Applied  Mechanics,  A.S.M.E.,  195B»  pp*  51-69" 

4.  L.  M.  Graves,  "The  theory  of  functions  of  real  variables".  Second 
Edition,  McGraw-Hill,  N.  Y.,  IG56,  Chapter  VIII,  Theorem  7,  pp.  150. 


-  28  - 


